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FOREWORD 


This  report  was  prepared  by  the  Department  of  Materials  Science 
and  Engineering,  Massachusetts  Institute  of  Technology.  The  report  is 
entitled,  "Electroslag  Remelting".  The  contract  was  administered  under 
the  technical  supervision  of  R.  H.  Frost  and  F.  C.  Quigley  of  the  Army 
Materials  and  Mechanics  Research  Center,  Watertown,  Massachusetts. 
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ABSTRACT 

This  interim  report  describes  research  conducted  at  the 
Massachusetts  Institute  of  Technology  on  macrosegregation  in  ingots 
produced  by  the  electroslag  remelting  (ESR)  process.  The  work  has  been 
on  producing  small  experimental  ingots  in  the  laboratory,  determining  the 
severity  of  segregation  in  those  ingots,  and  comparing  the  segregation 
with  calculations  using  macrosegregation  theory. 

Two  different  types  of  experimental  apparatus  are  described 
herein.  The  first,  employed  for  Al-4%Cu,  is  a small  scale  ESR  unit  which 
consists  of  a DC  power  source,  a water  cooled  mold,  consumable 
electrodes,  and  slag  layer  as  in  conventional  ESR  processing.  The  second 
type  of  apparatus  is  one  which  simulates  the  solidification  of  an  ESR 
ingot  but  does  not  employ  slag.  The  unit,  used  for  Sn-15%Pb,  comprises  a 
source  of  melt  droplets,  a cooled  mold,  and  a heat  source  to  simulate  the 
heat  input  of  the  ESR  process. 

Results  of  macrosegregation  in  a series  of  five  Al-4%Cu 
experimental  ingots  show  that  severe-localized  segregates,  as  sometimes 
found  in  large  commercial  ingots,  cannot  be  produced  using  this  alloy 
cast  in  the  small-scale  ESR  molds  (7.6  cm  in  diameter).  Surface  to 
center  compositional  variations  are  produced  enabling  a comparison  with 
calculations  using  macrosegregation  theory. 

In  order  to  study  a wider  range  of  segregation  problems 
encountered  in  ESR  ingots,  recent  efforts  are  directed  toward  solidifying 
Pb-15%Sn  alloy  in  the  simulated  ESR  apparatus.  Initial  results  obtained 
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are  very  encouraging.  Severe  macrosegregation  is  produced  in  this 
apparatus.  In  an  ingot  (8.2  cm  diameter),  the  composition  varies  from 
about  26%Pb  in  the  center  to  a minimum  of  6%Pb  slightly  beyond  midradius. 
Additionally,  there  are  "freckles"  at  approximately  midradius  and  an 
abnormally  high  fraction  of  eutectic  constituent  (similar  to 
"V-segregates")  in  the  center. 

The  analytical  work  is  based  upon  an  analysis  of  the  slow 
creeping  flow  of  liquid  between  dendrite  arms.  The  flow  is  driven  by 
(1)  gravity  induced  convection,  and  (2)  solidification  shrinkage. 

Agreement  between  calculations  of  macrosegregation  and  the  macrosegregation 
in  the  experimental  ingots  is  excellent  considering  the  large  number  of 
variables  involved.  The  theory  correctly  predicts  the  presence  of 
"freckles"  found  in  the  Sn-15%Pb  ingot. 
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I.  OVERVIEW  OF  RESEARCH  CONDUCTED  AND  RESULTS  OBTAINED 
A.  Background 

The  Casting  and  Solidification  Research  Group  at  the 

Massachusetts  Institute  of  Technology  has  studied  solidification  of  steel 

under  Army  sponsorship  for  about  sixteen  years. ^ Recent,  and 

current  work  is  on  the  processing  variables  which  affect  the 

solidification  structure  in  ESR  ingots.  A recent  report  described  work 

on  AISI  304  steel  made  in  specially  designed  experimental  ESR  molds. 

The  effects  of  power  input,  current  density,  slag  depth,  and  electrode 

diameter  on  electrode  melting  rate,  temperature  gradient  in  the  liquid, 

and  isotherm  velocities  were  studied  and  related  to  microstructure.  One 

innovation  resulting  from  that  work  is  that  by  electrode  and  slag  cooling 

the  investigators  were  able  to  achieve  increased  ratios  of  thermal 

gradient  to  solidification  rate  (G/R) . The  resulting  dendritic 

structures  were  much  finer  than  those  observed  in  exothermic  molds  at 

(12) 

comparable  distances  from  the  bottom  chill. 

Work  reported  herein  is  on  the  heat  and  fluid  flow  variables  which 

influence  macrosegregation  in  ESR  ingots.  Some  types  of 

macrosegregation,  including  "freckles”  and  surface  to  center  variations 

in  composition,  result  from  causes  that  can  be  described  quantitatively. 

The  M.I.T.  solidification  research  group  has  carried  out  numerous 

(13-19) 

analytical  and  experimental  studies  on  macrosegregation.  This 

interim  report  describes  work  which  applies  the  principles  of 
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macrosegregation  theory  to  the  solidification  of  ESR  ingots.  A computer 
analysis  is  developed  for  the  special  geometry  and  heat  flow  conditions 
found  in  experimental  ingots,  and  results  of  the  analysis  are  compared 
with  the  experimental  work  on  model  alloys,  Al-4%Cu  and  Sn-15%Pb. 

B . Experimental  Work 

Two  different  types  of  experimental  apparatus  were  employed  in 
this  work.  The  first  is  a small  ESR  unit  comprised  of  a DC  power  source, 
a water  cooled  mold,  a consumable  electrode  and  a slag  layer  as  in 
conventional  ESR,  Figure  1.  This  unit  was  used  for  experimental  work  on 
Al-4%  copper  alloy. 

The  second  apparatus  employed  is  one  which  simulates  solidification 
behavior  of  the  ESR  process  but  does  so  without  a consumable  electrode  or 
slag.  This  unit  comprises  a source  of  molten  metal  droplets,  a cooled 
mold,  and  a heat  source  to  simulate  ESR  heat  input.  Figure  2.  This  unit 
was  used  for  experimental  work  on  a Sn-15%Pb  alloy. 

In  both  of  the  apparati,  liquidus  and  solidus  shape  and  movement 
could  be  varied  significantly  and  could  be  made  to  closely  simulate  the 
behavior  of  actual  ESR  ingots.  For  example,  Figures  3 and  4 show  the 
shapes  of  isotherms  during  solidification  of  ingots  (one  of  Al-4%  Cu  and 
one  of  Sn-15%  Pb) . In  both  of  these  ingots,  significant  segregation  was 
found.  Figure  5 shows  the  segregation  in  the  Al-Cu  ingot  and  Figure  6 
the  segregation  in  the  Sn-15%  Pb  ingot.  The  segregation  in  the  Sn-15%  Pb 
ingot  was  significantly  severe  that  gross  segregates,  commonly  called 
"freckles",  formed,  as  shown  in  Figure  7. 
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C.  Analytical  Work 

An  underlying  concept  of  this  work  is  that  macr ©segregation,  such 

as  that  shown  in  Figures  5-7,  can  be  understood  and  quantitatively 

(13-19) 

analyzed  using  theories  developed  at  M.I.T.  over  the  last  decade. 

The  cause  of  the  segregation  is  believed  to  be  the  slow  creeping  flow  of 
liquid  between  dendrite  arms;  that  flow  is  driven  primarily  by 
(1)  gravity  induced  convection,  and  (2)  solidification  shrinkage. 

Basic  mathematics  describing  this  flow  was  developed  earlier  and 
is  summarized  herein.  In  this  work,  the  basic  equations  are  applied  to 
the  ESR  process  and  a computer  analysis  developed  to  solve  for  the  fluid 
flow  and  resulting  segregation  behavior  in  ESR  ingots.  As  examples  of 
the  computer  output.  Figures  8 and  9 show  the  calculated  interdendritic 
fluid  flow  in  the  two  ingots  referred  to  earlier.  With  flow  curves 
such  as  those  of  Figures  8 and  9,  calculations  of  macrosegregation  are 
readily  made.  Figure  10a  compares  the  experimental  results  of  the 
Al-4%Cu  ingot  with  theory,  and  Figure  10b  compares  the  results  of  the 
Sn-15%Pb  with  theory.  Agreement  is  excellent  considering  the  large 
number  of  variables  in  the  complex  segregation  process.  Detailed 
discussion  is  given  in  the  report  of  this  comparison  between  experiment 
and  theory,  as  well  as  of  the  choice  of  the  adjustable  parameter,  the 
’’permeability" , of  the  mushy  zone. 

In  addition  to  the  comparison  of  theory  with  experiment  described 
above,  it  should  be  noted  that  the  macrosegregation  theory  predicts  when 
"freckles"  such  as  those  shown  in  Figure  7 appear,  and  this  theory 
correctly  predicts  the  presence  of  the  freckles  in  that  ingot. 
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Figure  1:  Schematic  of  the  experimental  set-up  used  to  study 

macrosegregation  in  Al-4%  Cu  ESR  ingots. 
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Figure  3:  The  position  of  liquidus  and  solidus  isotherms  after 

7 minutes  in  ingot  1 (Al-4%  Cu) . 
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Figure  4:  The  shape  of  the  mushy  zone  in  the  Sn-15%  Pb  ingot  after 

14  minutes. 
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Figure 


Flow  lines  in  ingot  1 (Al-4.4%  Cu)  determined  using 
Y = 8 x 1CT8  cm2.  The  spacing  of  the  flow  lines  is 
approximately  proportional  to  the  inverse  of  the 
velocity  magnitude. 


HEIGHT 


16 


II.  INTRODUCTION 


A.  Related  Research 


The  ESR  process  is  one  of  the  most  important  new  processes 
developed  for  special  purpose  alloys.  The  main  advantages  are  the 
refining  that  can  be  obtained  by  melting  through  a slag  of  controlled 
composition,  and  the  special  control  over  solidification.  Such  control 
reduces  dendrite  arm  spacing,  microsegregation,  and  porosity,  hence 
leading  to  a sound  ingot. 

However,  the  recent  production  of  ESR  ingots  for  large  forgings 
has  shown  that  serious  macrosegregation  problems  including  centerline 


(21) 

segregation  and  channel-type  segregation  (freckles)  can  result. 

The  non-uniformity  of  properties  and  structure  can  deleteriously  affect 

the  mechanical  behavior  of  ingots  during  forging  or  rolling.  The 

development  of  the  continuous  ESR  process  is  also  limited  by  the  same 

(22) 

macrosegregation  problem. 

Like  other  solidification  processes,  the  shape  and  depth  of  the 
’’mushy1’  zone  piny  an  important  role  in  the  macrosegregation  found  in  ESR 
ingots.  It  has  been  found  that  the  operating  conditions  that  favor  a 
deeper  shape  of  mushy  zone  will  cause  more  macrosegregation  in  the  ESR 
ingots. xhe  reason  for  this,  according  to  the  macrosegregation 
theory^17, is  that  solute-rich  interdendritic  fluid  flow  can  be 
strongly  affected  by  solidification  shrinkage  and  gravity  in  a deep 


mushy  zone. 
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(10  25  26) 

Experimental  and  computer  studies^  5 * } have  shown  that  the 

melting  rate  of  the  electrode  is  the  most  important  parameter  in 
determining  the  depth  of  liquid  metal.  A higher  melting  rate  results  in 
a deeper  metal  pool  and  hence  a deeper  mushy  zone.  Other  operating 
parameters  such  as  composition  and  amount  of  slag,  thermal  conductivity 
of  metal,  current  and  voltage,  electrode  polarity  also  affect  the  depth 
of  liquid  pool  and  mushy  zone.  The  computer  has  also  been  used  to 
investigate  heat  flow  and  gap  formation  in  order  to  predict  conditions 
which  lead  to  break-outs  in  continuous  casting. 

Heat  flow  models  have  been  developed  to  relate  the  shape  of  mushy 

/fj/-  r\  q \ 

zone  to  the  operating  parameters  mentioned  above.  Up  to  now,  the 

only  quantitative  study  on  macrosegregation  in  ESR  ingots  was  carried 

out  by  Mitchell.  However,  his  approach  was  not  rigorous  since  no 

(31) 

macrosegregation  theory  was  considered.  Recently,  Thomas  et  al 

applied  a magnetic  field  shield  to  solidifying  ESR  ingots  and  reduced 

(21) 

macrosegregation.  Cooper'  y reported  on  the  so-called  "central  zone 

remelting  technique”  by  trepanning  and  remelting  the  badly  segregated 

central  zone  of  a conventionally  cast  big  ingot.  This  technique, 

however,  provides  no  direct  prevention  of  macrosegregation  formation  in 

the  final  core  solidified  as  an  ESR  ingot. 

As  the  solidification  research  group  at  M.I.T.  worked  on 
(13-19) 

macrosegregation  , other  researchers  reported  on  the  importance 

of  convection  in  the  mushy  zone  during  the  solidification  of  castings 

(24  32-34)  (34) 

and  ingots.  9 In  particular,  Copley  et  al  have  pointed  out 
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how  "density  inversions"  can  lead  to  the  formation  of  "freckles"  in 
nickel-base  super alloys. 

One  difficulty,  in  studying  the  solidification  of  ESR  ingots,  is 

(35) 

obtaining  thermal  data.  Mellberg  and  Sandberg v demonstrate  that  by 
measuring  secondary  dendrite  arm  spacing  and  "doping"  to  obtain  the 
shape  of  the  liquidus  isotherm,  it  is  possible  to  deduce  temperature 
distribution,  cooling  rates,  and  thermal  gradients  within  the  mushy  zone. 
Very  recently  Frost v has  applied  similar  technique  and  analyses  to 
describe  the  solidification  of  a series  of  AISI  4340  ESR  ingots  produced 
at  different  melt  rates. 

B . Permeability 

In  order  to  calculate  extent  of  macrosegregation  in  ESR  ingots, 
it  is  necessary  to  analyze  the  velocity  of  interdendritic  liquid  in  the 
mushy  zone.  DTArcyfs  law  is  used  to  describe  flow  through  porous 
media;  it  is 

v = - (VP  + piS}  (1) 

where  v is  the  interdendritic  liquid  velocity,  \i  is  the  viscosity,  g^ 
is  the  volume  fraction  liquid,  P is  pressure,  is  the  liquid  density, 
g is  the  acceleration  constant,  and  K is  the  specific  permeability  or 
simply  permeability.  Permeability,  as  defined  in  Equation  (1),  has 
dimensions  of  length  squared. 
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Physically,  the  permeability  can  be  pictured  in  two  ways,  both 
leading  to  the  same  result,  but  arrived  at  by  slightly  different 
physical  models.  In  one  case,  interdendritic  flow  is  pictured  as  taking 
place  in  long  tortuous  channels  through  the  porous  solid,  and  to  account 
for  the  fact  that  the  channels  are  neither  straight  nor  smooth,  a 
"tortuosity  factor",  T,  is  applied.  This  model  is  referred  to  by 
Mehrabian  et  alv  with  this  model  the  permeability,  K,  is 

2 

gT 

K = — (2) 

87TnT 

where  n is  the  number  of  flow  channels  per  unit  area  normal  to  flow,  and 
T is  the  tortuosity  factor.  For  a given  dendrite  arm  spacing,  if  n and 
T are  constants,  the  permeability  is  related  to  volume  fraction  liquid 
by 

K = Yg l (3) 

where 

y = (8TTnT^) 

(37) 

Piwonka  and  Flemings  report  values  of  permeability  measured  in 

Al-4%  Cu  alloy  which  obey  Equation  (3)  for  weight  fraction  liquid  up  to 

_9  2 

about  0.35  with  y - 6 x 10  cm  , Using  a somewhat  different 

( 38) 

experimental  technique,  Apelian  et  alv  measured  permeability  in 

porous  dendritic  networks  of  Al-4%  Si  alloys,  and  also  observed  that 

Equation  (3)  is  valid  for  volume  fraction  liquid  up  to  0.37  with 
-92  -92 

y ^ 9 x 10  cm  and  y - 3 x 10  cm  for  nongrain-refined  and 


20 


f37  38} 

grain-refined  structures,  respectively.  In  these  works k 9 
permeability  was  not  rigorously  determined  as  a function  of  dendrite  arm 
spacing.  Apelian  et  al  report  approximate  secondary  dendrite  arm 
spacings  of  280  microns  (nongrain-refined)  and  150  microns  (grain- 
refined)  measured  on  two  samples  only.  Piwonka  and  Flemings  do  not 
report  dendrite  arm  spacings,  but  it  appears  that  the  secondary  dendrite 
arm  spacing  is  about  140  microns.  The  results  of  these  measurements  of 
permeability  in  aluminum  alloys  are  shown  in  Figure  11.  For  g^  > 0.35, 
the  work  of  Piwonka  and  Flemings  show  that  the  slope  changes  to  a value 
much  greater  than  2 so  that  y in  Equation  (3)  is  not  constant.  This  can 
be  interpreted  to  indicate  that  the  number  of  channels  available  for 
flow  (n)  or  tortuosity  (t)  or  both  are  not  constant  when  g > 0.35. 

Li 

Also,  one  would  expect  in  the  range  where  y is  constant  that  permeability 

increase  with  increasing  dendrite  arm  spacing.  Figure  11  indicates  that 

this  is  not  necessarily  the  case,  although,  as  mentioned,  the  secondary 

arm  spacings  are  not  systematically  recorded  in  these  studies. 

A similar  model  for  permeability  is  to  imagine  flow  through 

channels  which  are  neither  straight  nor  smooth  with  areas  normal  to  flow 

which  are  highly  irregular.  In  this  case  the  channels  are  characterized 

(39) 

by  their  "hydraulic  radius",  R , which  is  shown  to  be 

H 


rh  ■ h,s 


(4) 


where  S is  the  surface  area  of  the  solid  phase  per  unit  volume  of  alloy 
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(i.e.,  per  unit  volume  of  "mush").  In  this  case. 


(5) 


and  a is  constant  in  the  range  of  4.2  - 5.  Equation  (5)  has  the 
advantage  that  in  order  to  calculate  permeability,  only  one  paramater 
of  the  structure,  S,  needs  to  be  characterized  whereas  Equation  (2) 
requires  knowledge  of  n and  T. 

Both  expressions.  Equations  (2)  and  (5),  can  be  used  to  predict, 
qualitatively,  the  effect  of  cooling  rate  during  solidification  on 
permeability.  In  ingots  with  high  cooling  rates,  dendrite  arms  are 
more  closely  spaced  than  in  ingots  slowly  cooled.  Since  n,  in 
Equation  (2),  increases  with  decreasing  dendritic  spacings,  then 
permeability  should  be  less  as  ingot  cooling  rate  increases.  Also,  the 
surface  area  of  the  dendritic  phase,  S,  would  increase  with  decreasing 
dendritic  spacing,  and  Equation  (5)  likewise  predicts  that  K should 
decrease . 

S treat  and  Weinberg^^  measured  permeability  in  partially 
solidified  Pb-20%  Sn  at  193°C.  In  their  experiments  volume  fraction 
liquid  was  held  constant  (g  = 0.19),  and  dendrite  arm  spacing  was 

Lj 

varied.  Using  average  primary  arm  spacings  (X),  and  assuming  that 

n = (A)“2  (6) 

they  tested  measured  values  of  K against  Equation  (2) . As  predicted  by 

_2 

Equations  (2)  and  (6) , permeability  varies  directly  with  X and  the 
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tortuosity  factor  is  4.6. 

Permeability  also  varied,  qualitatively,  as  expected  with 

secondary  arm  spacing,  but  no  simple  relationship  such  as  Equation  (6) 

could  be  readily  applied,  so  primary  spacings  were  considered  as  the 

more  useful  parameter  for  characterizing  the  structure However,  one 

could  argue  that  since  secondary  arm  spacing  relates  solely  to  cooling 

rate  during  solidification,  it  is  better  to  characterize  structures  by 

secondary  arm  spacings  even  though  a simple  relationship  between 

permeability  and  secondary  arm  spacing  apparently  does  not  exist.  The 

results  of  Streat  and  Weinberg are  given  in  Figure  12  in  terms  of 

secondary  arm  spacings.  For  the  largest  spacing  (83  microns)  examined, 

-9  2 

permeability  is  about  7 x 10  cm  ; extrapolating  results  to  arm 

_ & _ ~j  2 

spacings  of  140  microns,  permeability  is  about  10  - 10  cm  . At 

-10  2 

fraction  liquid  0.19,  permeability  in  aluminum  alloys  is  about  10  cm 
(Figure  11).  It  appears,  therefore,  that  the  permeability  (at  a given 
secondary  arm  spacing)  is  one  or  two  orders  of  magnitude  greater  in 
Pb-Sn  alloys  than  in  aluminum  alloys . 

To  sum  up,  studies  to  date  show  that  permeability  increases  with 
increasing  fraction  liquid  and  with  increasing  dendrite  arm  spacings 
but,  as  yet,  there  is  no  one  relationship  which  takes  into  account  both 
of  these  factors  quantitatively.  Because  of  these  uncertainties  in 
computing  permeabilities,  calculations  done  on  ingots,  reported  herein, 
use  permeability  as  given  by  Equation  (3) . 


LOG  PERMEABILITY,  cm 
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Figure  11:  Permeability  of  aluminum  alloys.  Numbers  on 

curves  are  approximate  secondary  arm  spacings. 


PERMEABILITY,  cm 


SECONDARY  DENDRITE  SPACING, /x 

Permeability  of  Pb-20%  Sn  alloy.  Experimental 
data  of  Streat  and  Weinberg^0) ; broken  line  is 
used  for  extrapolation. 


Figure  12: 
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III.  analysis;  of  macrosegregatiol 


A.  Theoretical  Model 


The  interdendritic  liquid  within  the  mushy  zone  is  acted  upon  by 
gravity;  consequently  if  the  density  of  that  liquid  varies  during 
solidification,  there  is  convection.  Ke  refer  to  this  as  "gravity- 
induced1'  convection.  Convection  of  interdendritic  liquid  also  occurs 
due  to  solidification  shrinkage  since  liquid  "feed  metal"  must  flow 
towards  regions  within  the  mushy  zone  where  the  solidifying  solid  has  a 
density  greater  than  the  local  interdendritic  liquid.  This  contribution 
to  the  convection  is  called  the  "solidification  induced"  convection. 

The  combined  effect  of  the  "gravity  induced"  and  the  "solidification 
induced"  convection  gives  the  total  convective  pattern  within  the  mushy 
zone  of  a solidifying  ingot.  If  solidification  time  is  long  (as  in  a 
large  steel  ingot),  "gravity  induced"  convection  predominates  because 
the  dendrites  in  the  solid-liquid  region  are  coarse  and  open  so  that 
solute-rich  interdendritic  liquid  flows  with  little  resistance  offered 
by  the  solid  phase.  Rather  extensive  flow  also  results  because  the 
time  available  for  flow  during  solidification  is  large.  In  a "chill" 
casting,  dendrites  are  fine  and  solidification  time  is  small;  hence, 
the  amount  of  flow  is  considerably  less  than  in  large  ingots  and  is 
primarily  "solidification  induced". 

To  describe  the  convection  of  interdendritic  liquid  within  the 
mushy  zone,  D'Arcy’s  Law  is  used,  so  that  Equation  (1)  applies. 
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The  interdendritic  liquid  velocity  v must  satisfy  continuity 
within  the  mushy  zone  which  means  that  at  all  times 

gs  + gL  = 1 (7) 

and 

If  (PsgS  + PLgL}  = "V'W  (8) 

where  g is  volume  fraction  of  the  solid.  Equation  (8)  is  the 
s 

continuity  equation  which  is  written  for  two  phases,  solid  and  liquid, 

with  the  assumption  that  the  solid  phase  is  stationary. 

Volume  fraction  liquid  varies  according  to  the  alloy  under 

consideration  and  the  manner  in  which  solute  partitions  between  the 

solid  and  liquid  phases  during  solidification,  the  thermal  gradient  VT 

and  cooling  rate  £,  and  the  amount  of  convection  within  the  mushy  zone. 

The  expression  which  takes  all  of  these  factors  into  account  has  been 

(13) 

derived  by  Flemings  and  Nereov  and  is  called  the  "local  solute 
redistribution  equation".  It  is 


Hi  + 


v*  VT 


(9) 


where  k is  the  equilibrium  partition  ratio,  3 is  the  solidification 

shrinkage  defined  as  (p  - pT)/p  and  C is  the  composition  of  the 

S h S Li 

interdendritic  liquid. 

Equations  (1)  and  (7  - 9)  are  combined  to  give  an  expression  for 
the  pressure  distribution  within  the  mushy  zone.  In  cylindrical 
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coordinates  (r  ,z)  , and  with  the  mushy  zone  moving  at  a constant  velocity, 
the  pressure  is  described  by 


+ A M + + B II 

3r2  + A3r+9z2  + "3z 


4-  C = 0 


(10) 


The  method  of  arriving  at  Equation  (10)  is  given  in  Appendix  A where  the 
coefficients  A,  B,  and  C,  are  also  defined.  Equation  (10),  along  with 
the  following  boundary  conditions,  constitute  the  pressure  distribution 
within  the  mushy  zone. 

By  symmetry,  there  is  no  radial  component  of  velocity  at  the 
centerline  so  the  pressure  gradient  in  the  r-direction  is  zero;  therefore 


§7  (z *0)  = o . (u) 

Also,  by  virtue  of  the  mold  wall  there  is  no  radial  velocity  component 
at  the  ingot  wall  (r=R) ; so 


f (z,R)  = 0 


(12) 


Within  the  bulk  liquid,  we  assume  no  flow  so  that  at  the  liquidus 
isotherm  the  pressure  is  that  of  a motionless  liquid;  i.e., 


P0  + 


PL0gh 


P (liquidus) 


(13) 
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where  PQ  is  the  pressure  at  the  top  of  the  liquid  pool  (approximately 

1 atm.),  PTn  is  the  density  of  the  liquid,  g is  the  gravitational 
LU 

constant,  and  h is  the  height  of  the  liquid  pool  which  varies  with  radius 
unless  solidification  is  unidirectional. 

At  the  solidus  isotherm,  liquid  feeds  the  contraction  of  the 
solidifying  eutectic  so  that 

v (solidus)  = - ( PsE  — ^-E  ) U (14) 

PLE  E 

U is  the  solidus  (eutectic)  isotherm;  in  general,  its  components  vary 
E 

along  the  solidus.  Solidified  eutectic  density  is  and  the  density 

of  the  eutectic  liquid  is  PT„. 

Lit, 

To  solve  Equation  (10)  requires  that  the  temperature  distribution 
within  the  mushy  zone  be  known.  In  this  work,  temperature  is  determined 
experimentally.  With  temperature,  the  coefficients  in  Equation  (10)  can 
be  determined  because  the  local  composition  of  the  interdendri tic  liquid 
is  only  a function  of  the  temperature  as  given  by  the  phase  diagram. 

For  the  model  system  employed,  Al-4%  Cu  and  Sn-15%  Pb,  the  phase  diagrams 
used  are  given  in  Figures  13a  and  14a.  It  is  also  necessary  that  the 
local  values  of  liquid  and  solid  densities  be  known  throughout  the  mushy 
zone;  the  densities  are  given  in  Figures  13b  and  14b  for  the  two  model 
alloys.  Finally,  the  local  volume  fraction  of  liquid,  must  also  be 

known  to  calculate  the  coefficients  in  Equation  (10);  however,  their 
values  are  not  known  until  Equation  (10)  is  first  solved  and  then 
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Equation  (9)  is  satisfied.  Therefore,  a numerical  technique  is  employed 

to  arrive  at  the  correct  values  of  g^. 

The  sequence  of  calculations  directed  by  the  computer  code  is 

outlined  in  Figure  15.  After  the  input  has  been  assimilated,  the  pressure 

distribution  is  found  according  to  f inite-dif ference  approximations  of 

Equations  (10)  - (14).  At  this  point,  reasonable  approximations  of  gL 

(47) 

are  used  as  input;  we  have  used  the  Scheil  equation  in  the  form 
which  gives  g as  a function  of  temperature.  The  velocity  field  within 

Li 

the  mushy  zone  is  then  evaluated  according  to  Equation  (1)  which  leads 
directly  to  a solution  of  Equation  (9)  giving  "updated"  values  of  g^. 
Calculations  are  repeated  by  replacing  the  initial  values  of  g (estimated 
with  the  Scheil  equation)  with  "updated"  values.  The  process  is 
continued  until  there  is  no  change  in  values  of  gL  in  the  mushy  zone. 

Finally,  macrosegregation  in  an  ingot  is  given  in  terms  of  the 


In  Equation  (15),  g is  the  volume  fraction  of  eutectic  and  C is  the 

E ti 

eutectic  composition.  The  integration  in  the  numerator  can  be  carried 
out  by  picturing  a fixed  unit  volume  in  which  C continuously  changes  in 
a known  manner  as  the  alloy  solidifies.  For  an  ESR  ingot,  when  the 
isotherms  move  at  a steady  velocity,  the  integration  can  be  carried  out 
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by  integrating  from  the  liquidus  (g  = 0)  down  to  the  solidus  (g  = 1-g  ) 
at  a given  radius.  By  doing  this  for  different  radial  positions  within 
the  ingot,  versus  radius  is  determined  which  can  be  plotted  to  give 
the  pattern  of  macrosegregation. 

To  sum  up,  the  analytical  model  (given  a temperature  distribution) 
can  be  used  to  predict: 

(1)  the  pressure  distribution  within  the  mushy  zone; 

(2)  the  velocity  of  interdendritic  liquid  flow  within  the  mushy 
zone; 

(3)  the  distribution  of  volume  fraction  liquid  within  the  mushy 
zone;  and 

(4)  the  local  average  composition  after  solidification  is  complete. 
Accordingly,  from  step  (4)  the  extent  and  pattern  of  macrosegregation  is 
calculated  which  shows  a continuous  variation  in  composition  from  surface 
to  center.  Appendix  B discusses  in  more  detail  the  finite  difference 
equations  and  other  features  of  the  computer  code  used  to  analyze 
macrosegregation . 

B.  Flow  Instability  and  Localized  Segregates 

Up  to  this  point  solidification  is  pictured  as  progressively  taking 
place  as  the  temperature  within  the  mushy  zone  decreases.  Interdendritic 
liquid  flows  towards  regions  of  increasing  amounts  of  solid;  i.e.,  flow 
is  from  hotter  to  colder  regions  in  the  mushy  zone  and  is  said  to  be 
’’stable".  However,  as  discussed  by  Mehrabian  et  al^^ 


, under  certain 
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conditions  flow  may  become  "unstable"  with  consequent  formation  of  highly 
segregated  channels  called  "freckles"  when  they  occur  in  ESR  ingots. 

Figure  16  shows  flow  channels  under  conditions  of  "intermediate" 
flow,  and  "unstable"  flow.  These  flows  occur  whenever  vVT/e  < 0;  that  is, 
whenever  flow  is  from  cooler  to  hotter  regions.  This  can  most  easily  be 
seen  by  first  multiplying  Equation  (10)  by  (dC  /dT)  to  obtain  the  local 

Lt 

fraction  of  liquid  as  a function  of  temperature  which  is 


(1  + 


v-VT, 

e ' 


(16) 


The  quantity  inside  the  brackets  is  always  negative,  regardless  of  the 
value  of  k.  Now,  imagine  a perturbation  such  that  one  of  the  channels  in 
the  volume  element  sketched  in  Figure  16  becomes  slightly  larger  than  the 

surrounding  channels.  Resistance  to  flow  through  this  channel  is 

-)■ 

decreased,  and  flow  (absolute  value  of  v)  is  increased,  whatever  its 
direction.  If  in  addition,  however,  flow  is  in  a direction  of  increasing 
temperature  (v*VT/e  < -1)  then  Bg  /9T  is  negative,  i.e.,  as  temperature 

Lj 

decreases,  the  local  volume  fraction  liquid  increases  and  so  there  is  a 
driving  force  for  the  formation  of  channels;  local  melting  occurs. 

Consider  the  effect  of  a minor  perturbation  in  the  size  of  a flow 
channel  when  some  flow  occurs  in  the  direction  of  isotherm  movement. 
Specifically,  consider  the  case  when  0 < v'VT/e  < -1.  From  Equation  (16), 
the  greater  the  flow,  the  less  is  the  rate  of  solidification  in  a given 
channel.  As  sketched  in  Figure  16,  therefore,  the  size  of  the  perturbed 
channel  does  not  decrease  as  rapidly  as  that  of  its  neighbors  during 
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solidification.  However,  as  long  as  vVT/e  remains  within  the  limits  given 
above,  3g^/3T  remains  negative,  and  so  no  remelting  occurs.  If,  however, 
the  condition  vVT/e  < -1  is  reached  in  the  perturbed  channel,  remelting 
does  occur,  and  the  size  of  the  perturbed  channel  actually  grows  during 
cooling  of  the  volume  element.  Growth  can  occur  until  it  becomes  optically 
visible  as  a "freckle”. 
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(a) 


(b) 

Phase  diagram  and  density  used  for  calculations  of 
macrosegregation  in  Al-4%  Cu  ESR  ingots.  (a)  Phase 
diagram.  Ref.  41;  (b)  density  of  solid  and  liquid 

during  solidification.  Ref.  17. 


Figure  13: 
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WEIGHT  PERCENT  LEAD 


(a) 


(b) 

Figure  14:  Phase  diagram  and  density  used  for  calculations  of 

macrosegregation  in  Sn-20%  Pb  ESR  ingots.  (a)  Phase 
diagram.  Ref.  42;  (b)  density  of  solid  and  liquid 

during  solidification.  Refs.  43-46. 
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STABLE  II 

Figure  16:  Schematic  illustrat: 
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IV.  APPARATUS  AND  EXPERIMENTAL  PROCEDURE 


A.  General 

Two  different  types  of  experimental  apparatus  were  employed  in  this 
work,  and  these  are  each  described  below.  The  first  (employed  for  Al-4% 

Cu  alloy)  is  a small  scale  ESR  unit.  This  unit  comprises  a DC  power 
source,  a water  cooled  mold,  a consumable  electrode,  and  slag  layer  as  in 
conventional  ESR,  Figure  1. 

The  second  type  of  apparatus  employed  is  one  which  simulates  the 
solidification  behavior  of  the  ESR  process,  but  does  not  employ  slag. 

This  unit,  used  for  Sn-15%  Pb  alloy,  comprises  a source  of  molten  alloy 
droplets,  a cooled  mold,  and  a heat  source  to  simulate  the  heat  input  of 
the  ESR  process.  Figure  2. 

B.  Small  Scale  ESR  Unit 

1 . Apparatus 

A sketch  of  the  apparatus  is  shown  in  Figure  1 and,  with  the 

(12) 

exception  of  the  mold,  is  the  same  previously  used  by  Basaran  et  al 

The  power  supply  is  a D.C.  arc-welder  capable  of  providing  up  to  1600  amp. 

and  40  volts  for  a total  power  of  55  KW  (Miller  Electric,  SR  1000  BI) . 

The  electrode  mount  is  connected  to  a feed  screw  about  6 feet  long  which 
is  driven  by  one  of  two  gear  reduction  boxes  in  series  with  an  electric 
motor.  The  driving  speed  of  the  electrode  is  controlled  by  a Minarik 
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controller  (Model  SH-63);  in  this  way  the  electrode  can  be  operated  with  a 
linear  speed  in  the  range  of  0.2  - 14  cm/min.  The  electrodes  used  were 
2024  rods  (Al-4%  Cu)  1 inch  in  diameter  by  6 feet  long. 

In  the  initial  stages  of  the  work,  copper  molds  were  designed  and 
fabricated  to  provide  for  the  flow  of  cooling  water  through  the  annulus 
between  the  closed  jacket  and  the  mold  exterior  (Figure  17) . The  inside 
dimensions  of  the  mold  were  3 inches  in  diameter  by  8 inches  high.  A 
thin  wash  of  aluminum  oxide  powder  was  applied  to  the  mold  wall  in  an 
attempt  to  eliminate  corrosion  of  the  mold  by  the  slag  (45%  LiCl  - 55% 

KC1) . Several  preliminary  ingots  were  made  in  order  to  establish  workable 
electrode  power  settings,  slag  volumes,  and  start-up  techniques.  Mold 
corrosion  was  a persistent  problem  particularly  since  the  corrosion  was 
sometimes  localized  making  it  necessary  to  fabricate  new  molds.  Other 
mold  washes  were  tried,  such  as  a graphite  wash  and  a wash  composed  of 
graphite  and  zirconia.  The  latter  reduced  the  severity  of  the  slag 
attack  but  molds  were  still  attacked  perchance  the  coating  spalled  during 
a run. 

A second  mold  was  designed  mainly  to  incorporate  the  feature  of 
being  able  to  merely  replace  the  mold  wall  rather  than  completely  fabricate 
an  entire  mold.  This  design  also  made  it  easy  to  insert  up  to  5 
thermocouples  into  the  ingots  during  a run.  This  design  is  shown  in 
Figure  18.  Cooling  water  was  sealed  from  the  mold  interior  at  the 
thermocouple  entries  with  rubber  rings.  The  replaceable  copper  mold  was 
3 inches  in  diameter  by  14  inches  in  height.  An  aluminum  bottom  hearth 
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was  used  for  bottom  welding  at  start-up  and  to  prevent  slag  corrosion  of 
the  copper  bottom  chill. 

2.  Experimental  Procedure 

About  250  c.c.  of  liquid  KC1  - LiCl  slag  of  eutectic  composition 
(45  wt%  LiCl  and  55  wt%  KC1)  was  poured  into  the  mold  at  start-up. 
Immediately  after  pouring  the  slag,  the  arc-welder  was  turned  on  to 
initiate  melting  with  a power  of  10  KW.  The  cooling  water  was  turned  on, 
and  then  the  power  was  lowered  to  the  working  value  between  2 and  4 KW. 
The  amperage  was  kept  constant  by  adjusting  the  driving  speed  of  the 
electrode,  and  the  electrode  position  was  recorded  during  a run.  When 
thermal  data  were  desired,  five  chromel-alumel  thermocouples  were  used. 
They  were  inserted  into  the  mold  and  pushed  to  predetermined  positions 
(shown  as  X!s  in  Figure  18)  immediately  after  the  electrode  passed  these 
positions,  and  their  output  was  recorded.  Cooling  curves  obtained  were 
used  to  determine  the  shape  of  mushy  zone,  temperature  distribution, 
cooling  rates  and  solidification  rates.  Thermal  data  for  two  ingots 
(Numbers  1 and  2)  were  obtained  with  the  second  mold  design  (Figure  18). 

Three  ingots,  made  in  the  first  mold  design  (Figure  17)  were 
doped  with  about  five  grams  of  Al-50%  Cu  to  reveal  the  liquid  pool  shape. 
The  ingots  were  cut  into  sections  in  order  to  obtain  analyses  by  X-ray 
fluorescence  to  detect  macrosegregation  across  the  ingots.  The  ingots 
were  also  macroetched  with  Keller’s  etch  to  reveal  the  doping  effect, 
and  samples  for  microstructures  were  removed. 
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C.  Simulated  ESR  Apparatus 

1 . Apparatus 

A sketch  of  the  apparatus  is  shown  in  Figure  2.  The  stainless 
steel  mold  is  3-1/4  in.  in  diameter  and  13  in.  long.  The  metal  pool 
inside  the  mold  is  heated  with  six,  3 in.  long,  resistance  heaters 
(G  3A133  Wattlow  Electric  Manufacturing  Company)  connected  in  parallel. 
These  heaters  were  positioned  inside  holes  drilled  into  a 3 in.  long  by 
1-1/4  in.  diameter  stainless  steel  bar.  They  were  connected  in  series  to 
a 240-volt  variare  used  to  adjust  and  control  the  power  input. 

Cooling  water  runs  through  a movable  cooling  jacket  surrounding 
the  mold.  Both  the  resistance  heaters  and  the  cooler  are  fixed  to  the 
same  system  used  for  driving  the  electrodes  in  making  the  Al-4%  Cu 
ingots.  Thermal  measurements  are  made  with  three  chromel-chromel 
thermocouples  located  inside  three  vertical  stainless  steel  tubes.  The 
tubes  are  fixed,  but  the  position  of  the  thermocouples  is  varied  during 
a run  by  sliding  them  up  and  down  inside  the  tubes. 

In  this  analog  set-up,  the  flow  of  liquid  Sn-15%  Pb  alloy  from  the 
top  stainless  steel  container  is  controlled  by  an  adjustable  valve. 
Heating  of  the  melt  in  the  top  container  is  done  with  two  8 in.  diameter 
by  1.5  in.  wide  band  heaters  (7533EX,  Wattlow  Electric  Manufacturing 
Company) . These  heaters  are  controlled  by  a thermocouple  hooked  up  with 
a Eurotherm  temperature  controller  (917/FCP/K/600/115/230V/A/C2) . A 
stirrer  was  used  to  insure  uniform  concentration  distribution  in  the 
bath  of  the  top  container. 
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2.  Experimental  Procedure 

Pure  tin  (99.9%)  and  pure  lead  (99.9%)  of  known  weights  were 
melted  together  in  a crucible  furnace.  The  liquid  alloy  was  well  stirred 
and  about  one-fourth  of  the  charge  was  poured  into  the  stainless  steel 
mold  until  the  liquid  level  rose  to  almost  the  top  of  the  resistance 
heaters  (about  3.5  in.  from  the  bottom  of  the  mold).  The  remaining  alloy 
was  then  poured  into  the  top  container.  Both  the  band  heaters  and  the 
resistance  heaters  were  then  turned  on.  The  cooling  water  was  also 
turned  on  at  the  same  time. 

With  the  resistance  heaters  and  the  cooling  jacket  fixed,  the 
initial  position  and  the  shape  of  the  mushy  zone  were  determined  by  moving 
the  three  thermocouples  up  and  down  and  locating  the  position  of  the 
liquidus  and  solidus  temperatures  of  the  alloy.  Power  input  to  the 
heaters  and  cooling  jacket  position  were  adjusted  until  the  desired 
position  and  shape  of  mushy  zone  were  obtained. 

Then,  the  resistance  heaters  and  the  cooler  were  moved  upwards  at 
the  predetermined  speed  and  the  valve  of  the  top  container  was  adjusted 
so  that  the  liquid  level  inside  the  mold  rose  at  the  same  speed.  As 
solidification  progressed  the  three  thermocouples  were  individually 
moved  up  and  down  so  that  a trace  of  solidus  and  liquidus  temperature 
could  be  made  along  with  the  temperature  distribution  in  the  mushy  zone. 
The  position  and  time  of  each  thermocouple  reading  were  recorded.  These 
thermal  measurements  were  used  to  determine  the  position  and  the  shape 
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of  the  mushy  zone,  the  temperature  distribution,  and  the  cooling  rates 
during  solidification. 

The  ingot  cast  was  cut  into  sections  for  microstructural  study  and 
analysis  by  X-ray  fluorescence  to  determine  degree  of  macrosegregation. 

D.  Chemical  Analysis 

Chemical  analysis  of  macrosegregation  in  the  ingots  was  by  X-ray 
fluorescence.  A General  Electric  X-ray  diffraction  unit  (Model  XRD3 
Type  1)  was  used  with  a Mo  tube.  The  primary  white  radiation  from  the 
tube  fluoresced  the  samples  on  an  area  of  0.32  cm.  diameter.  This  area 
covered  many  dendrite  arms  (secondary  dendrite  arm  spacing  is  about  40y) 
and,  therefore,  the  compositions  measured  were  local  average  compositions. 

The  secondary  radiation  from  the  sample  was  received  by  a Si  (Li) 
X-ray  detector  (ORTEC  Model  78915-10225)  and  the  intensity  of  the 
characteristic  line  (K^  for  Cu  and  for  Pb)  was  compared  with  a 
standard  intensity  versus  composition  curve  to  determine  the  composition. 

The  standards  were  prepared  from  rapidly  cooled  thin  sections  of 
known  compositions.  The  Al-Cu  standards  used  were  those  prepared  by 
Nereo^^.  The  Pb-Sn  standards  were  prepared  by  melting  pure  lead  and 
pure  tin  together  in  a graphite  crucible  to  form  the  liquid  alloy  of 
desired  composition.  The  liquid  alloy  was  heated  above  500°C,  well 
stirred  and  skimmed.  About  5 grams  of  the  liquid  alloy  was  quickly 
removed  from  the  crucible  and  dropped  1.5  feet  onto  a copper  chill  (1.5 
inches  x 6 inches  x 10  inches).  The  descending  drop  of  liquid  alloy  hit 
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the  copper  and  solidified  very  rapidly  as  a splat.  A significant  portion 
of  the  splat  formed  an  area  which  was  similar  to  a thin  disk  (about  0.8  mm. 
thick  and  2 in.  in  diameter)  with  a very  smooth  and  flat  bottom  surface. 
From  the  disk,  a thin  plate  was  then  cut  into  a rectangular  plate  (0.8  inch 
x 1.2  inches)  and  polished  to  600  grit  metallographic  paper  for  X-ray 
fluorescence.  The  remainder  of  the  disk  was  analyzed  by  wet  chemical 
analysis  to  determine  the  composition. 
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Figure  17:  ESR  mold  design  with  brazed  cooling-water  jacket. 
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Figure  18:  ESR  mold  design  with  replaceable  mold  wall.  Numbers 

indicate  thermocouples. 
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V,  EXPERIMENTAL  RESULTS  OF  MACROSEGREGATION 


A.  Aluminum-4%  Copper  ESR  Ingots 

Results  of  five  ingots  are  given  which  encompass  an  ingot  with  the 
deepest  isotherm  shapes  obtained  (No.  1),  an  ingot  with  flat  isotherms 
(No.  2;  i.e.,  unidirectional)  and  three  ingots  with  isotherms  of  shapes 
intermediate  between  the  deepest  and  the  flat  isotherms  (Nos.  3,  4,  and 
5).  The  isotherms  for  ingots  1 and  2 were  obtained  from  thermal 
measurements;  for  ingots  3,  4,  and  5 the  shapes  of  the  isotherms  were 
determined  by  doping. 

The  thermal  data  and  measured  compositions  (to  show  degree  of 
macrosegregation)  in  these  five  ingots  follow  in  this  section. 

Additional  discussion  of  the  experimental  results  is  given  later  in  this 
report  in  conjunction  with  calculations  resulting  from  the  analytical 
study  of  macrosegregation. 

1.  Ingot  No.  1 

The  cooling  curves  obtained  from  the  thermocouples  as  positioned 
in  Figure  18  are  given  in  Figure  19.  Thermocouples  1,  2,  and  3 were 
located  along  the  ingot  axis  respectively,  from  the  bottom  chill.  At  a 
given  temperature,  the  time  lag  between  thermocouples  1 and  2 is  somewhat 
greater  than  between  2 and  3 indicating  that  the  solidification  rate 
increased  slightly  as  solidification  proceeds  from  the  bottom.  As 
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expected,  thermocouples  4 and  5 cooled  sooner  than  No.  3 because  they  were 
located  closer  to  the  ingot  wall.  The  interesting  fact  regarding  these 
cooling'  curves,  however,  is  that  cooling  rate  over  the  solidification 
range  is  independent  of  the  distance  from  the  bottom  chill. 

Figure  20  shows  plots  of  the  position  of  the  liquidus  isotherm  (Z^) 
and  eutectic  isotherm  (Z  ) at  the  center  of  the  ingot.  After  steady 

iii 

state  is  achieved,  the  isotherms  move  with  a vertical  speed  of  0.053  cm/s. 
Figure  21  shows  that  the  isotherm  speed  in  the  axial  direction,  at  7.6  cm. 
up  from  the  bottom  chill,  is  independent  of  radius  since  the  time  lag 
between  the  liquidus  and  solidus  are  equal  for  the  three  thermocouples. 

Figures  19  - 21  can  be  used  to  determine  the  shape  of  the  mushy 
zone.  This  is  shown  in  Figure  3 after  7 minutes  of  ingot  solidification 
has  elapsed.  As  a consequence  of  cooling  rate  not  varying  with  radius, 
the  isotherms  are  parallel.  Finally,  temperature  gradient  in  the  mushy 
zone  is  shown  in  Figure  22.  Since  isotherms  are  parallel,  this 
temperature  gradient  is  independent  of  radius.  Figures  19  - 22  provide 
the  thermal  data  necessary  as  input  to  the  computer  simulation  which 
calculates  macrosegregation  in  Ingot  1. 

Figure  5 shows  the  degree  of  macrosegregation  obtained  in  Ingot  1 
as  measured  by  X-ray  fluorescence.  The  overall  analysis  of  the  ingot  is 
4.4%  Cu;  as  a result  of  solidification,  there  is  positive  segregation  at 
the  surface  (about  4.6%  Cu)  and  negative  segregation  at  the  center 
(about  4.25%  Cu) . Macroetching  showed  no  evidence  of  localized 
segregates,  such  as  "freckles"  or  "V"-segregates  along  the  center  line, 
so  segregation  is  continuous  from  center  to  surface.  Since  commercial 
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2024  aluminum  alloy  containing  titanium  was  used  as  electrode,  the  ingot 
structure  is  equiaxed  and  fine-grained;  this  is  true  of  all  five 
aluminum  ingots. 

Because  the  structure  is  equiaxed,  measurements  of  the  average 
dendrite  cell  size,  d,  were  made  rather  than  the  often  reported 
measurements  of  secondary  dendrite  arm  spacing.  Here  we  define  the 
dendrite-cell  size  as  the  distance  between  the  centers  of  adjacent  primary 
phases  separated  by  eutectic  constituent.  The  average  value  of  the 
dendrite-cell  size  is  obtained  by  drawing  random  lines  over 
photomicrographs  and  counting  the  numer  of  intercepts  between  the  lines 
and  the  eutectic  constituent.  About  200  intercepts  for  each  sample  were 
counted  in  photomicrographs  with  a magnification  of  128X.  Dendrite  cell 
size,  d,  in  Ingot  No.  1 is  shown  in  Figure  23.  No  significant  variation 
in  d is  observed  which,  of  course,  is  predictable  by  the  fact  that 
cooling  rate  during  solidification  is  also  constant  across  the  ingot. 

2 . Ingot  No.  2 

This  ingot  was  cast  and  examined  in  the  same  manner  as  Ingot  No.  1 
with  the  exception  that  an  unusually  thick  (about  0.3  cm.)  coating  of 
mold  wash  (graphite  and  powdered  zirconia)  was  applied  to  the  inside  mold 
wall.  As  a result,  solidification  was  unidirectional. 

Figures  24  and  25  show  cooling  curves  and  isotherm  positions, 
respectively.  At  7.6  cm.  from  the  bottom,  temperature  is  independent  of 
radius  (Figure  24)  so  isotherms  must  be  flat.  Cooling  rate  during 
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solidification  is  about  equal  to  that  observed  in  Ingot  1,  but  the  vertical 
temperature  gradients  are  somewhat  greater  (Figure  25) . Figure  26  shows 
the  mushy  zone  after  7 minutes  have  elapsed,  and  Figure  23  shows  the 
average  dendrite  cell  size  which  is  equal  to  that  observed  in  Ingot  No.  1. 

Since  solidification  is  unidirectional,  no  macrosegregation  is 
detected  in  this  ingot;  nor  is  there  any  evidence  of  localized  segregates 
found  in  an  etched  macrosection. 

3.  Ingots  3,  4,  5 

Thermal  data  for  these  ingots  were  not  measured,  but  with  isotherm 
shapes  obtained  by  doping  and  cooling  rates  calculated  using  measurements 
of  dendrite  cell  sizes,  the  thermal  history  can  be  constructed. 

Figure  27  shows  the  etched  macrostructure  of  Ingot  No.  4.  The 
important  feature  seen  is  that  doping  outlines  the  shape  of  an  isotherm. 
Since  doping  was  carried  out  by  suddenly  increasing  the  composition  of 
copper  in  the  bulk  liquid,  then  a melting  back  of  the  dendrite  tips 
protruding  into  the  liquid  pool  must  have  occurred.  Therefore,  the 
extent  of  this  remelting  must  be  known  in  order  to  predict  the  precise 
isotherm  indicated  by  the  doping  in  Figure  27.  However,  since  thermal 
data  for  Ingots  1 and  2 indicate  that  isotherms  are  parallel  and  Ingots  1 
and  2 represent  the  extremes  in  isotherm  shapes,  we  assume  that  isotherms 
in  Ingots  3,  4,  and  5 are  also  parallel.  As  such,  exactly  which 
isotherm  is  brought  out  by  doping  is  not  an  important  consideration,  but 
it  is  important  to  detect  the  shape  of  the  isotherms. 
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Another  indication  that  isotherms  are  parallel  in  these  ingots  is 
given  in  Figure  28;  dendrite  cell  size  in  these  ingots  is  uniform  across 
each'  of  the  ingots  which  is  consistent  with  the  existence  of  parallel 
isotherms. 

Figure  28  is  used  to  deduce  the  local  solidification  times  in  Ingots 

3,  4,  and  5.  For  aluminum  alloys,  the  secondary  arm  spacing,  d is  related 

(49) 

to  solidification  time,  0,  by  : 

d = 7.5  X 10-4e0-39  (17) 

with  d in  cm.  and  0 in  seconds.  Dendrite  cell  size,  d,  as  defined  herein 
would  only  equal  d when,  perchance,  the  plane  of  polish  includes  a row  of 
aligned  secondary  dendrite  arms  and  is  parallel  to  the  center  of  the  arms. 
In  general,  therefore,  for  random  planes  d is  greater  than  d,  but  the 
ratio  d/d  is  a constant  fraction.  As  such, 

d = A0°* 39  (18) 

in  which  A is  constant  and  greater  than  the  proportionality  constant  in 
Equation  (17) . 

Measured  local  solidification  time  of  Ingot  No.  1 is  95  s (average 

obtained  from  thermocouples  3,  4,  and  5,  in  Figure  19);  average  dendrite 

-4 

cell  size  is  61.5  x 10  cm.  (Figure  23).  Similarly,  for  Ingot  2 

-4 

(Figures  23  and  24),  0 is  78  s and  d is  61  x 10  cm.  By  calculating  A 
for  Ingots  1 and  2,  and  using  the  average  value,  we  arrive  at 

-3  0.39 

= 1.08  x 10  °0 


d 


(19) 
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Values  of  d obtained  from  Figure  28  are  then  substituted  into  Equation  (19) 
to  determine  the  local  solidification  times  in  Ingots  3,  4,  and  5. 

Results*  are  89,  89,  and  54  s,  respectively. 

Degree  of  macrosegregation  for  Ingot  3 is  shown  in  Figure  29.  The 
extent  of  segregation  is  considerably  less  than  that  detected  in  Ingot  1 
which  has  a deeper  isotherm  shape.  The  isotherms  are  more  nearly  flat  in 
Ingots  4 and  5 (Figure  30) , enough  so  that  no  segregation  is  detected  in 
these  ingots  (Figure  31). 

B.  Tin-15%  Lead  ESR  Analog  Ingot 

Results  of  macrosegregation  in  the  series  of  the  five  Al-4%  Cu 
experimental  ingots  show  that  severe-localized  segregates,  as  sometimes 
found  in  large  commercial  ESR  ingots,  cannot  be  produced  using  this  alloy 
cast  by  ESR  in  small  laboratory  scale  (3  inches  in  diameter)  molds. 

Although  surface  to  center  compositional  variations  were  produced,  these 
variations  are  rather  modest  and  no  severe-localized  segregates  were 
detected.  It  was  decided,  therefore,  to  design  laboratory  experiments 
which  could  be  used  to  study  a wider  range  of  segregation  problems 
encountered  in  ESR  ingots.  Accordingly,  experimental  efforts  were 
directed  towards  solidifying  Pb-15%  Sn  alloy  in  the  simulated  ESR 
apparatus  (Figure  2).  Results  obtained  from  our  initial  ingot  are  very 
encouraging  and  experimental  emphasis  will  continue  on  the  analog  apparatus 
and  the  Sn-Pb  alloy  selected  as  the  model  alloy. 
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Figure  32  shows  thermal  data  obtained  from  the  experiment  made  in 
the  ESR- analog  apparatus.  The  motion  of  the  liquidus  and  non- equilibrium 
solidus  is  plotted  for  three  radial  positions:  at  the  center,  slightly 

beyond  mid-radius,  and  almost  at  full-radius.  After  about  5 minutes  all 
curves  are  linear  and  parallel  which  indicates  steady  state  is  achieved 
above  9 cm.  from  the  ingot  base.  But  the  distance  between  isotherms 
varies  with  radius  so  the  mushy  zone  is  not  comprised  of  parallel 
isotherms  as  in  the  experimental  Al-4%  Cu  ESR  ingots.  Rather,  as  shown 
in  Figure  4,  the  mushy  zone  is  widest  at  ingot  center  and  decreases  as 
ingot  radius  is  approached.  Predictably,  the  microstructure  shows  finer 
dendrites  near  ingot  radius,  and  they  coarsen  as  the  ingot  center  is 
approached  (Figure  33) . 

In  this  ingot,  macrosegregation  is  severe.  Figure  6 shows  that  the 
composition  varies  from  a high  of  about  26%  Pb  in  the  center  to  a minimum 
of  6%  Pb  slightly  beyond  mid-radius.  The  curve  showing  the  overall 
segregation  can  only  be  drawn  approximately  since  there  are  localized- 
segregates  ("freckles”)  in  this  ingot.  Figure  7 shows  "freckles”  at 
approximately  mid-radius  and  an  abnormally  high  fraction  of  eutectic 
constituent  in  the  central  zone.  From  an  experimental  viewpoint, 
therefore,  the  analog  ESR  experiment  provides  a method  of  producing,  in 
the  laboratory,  the  severe  segregation  sometimes  found  in  large  ESR 
ingots.  As  discussed  later,  our  analytical  studies  confirm  this 
conclusion  and  so  we  plan  to  conduct  future  experiments  using  this 
apparatus.  Currently,  the  apparatus  is  being  used  by  one  of  the  authors 
(S.K.)  to  obtain  experimental  results  for  his  Ph.D.  thesis  research  work. 


Figure  20:  Liquidus  and  solidus  isotherms  along  axis  in  ingot 

(Al-4%  Cu). 
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Figure  21:  Liquidus  and  solidus  isotherms  with  respect  to  radius  at 

7.6  cm  from  the  bottom  in  ingot  1 (Al-4%  Cu) . 
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Figure  22: 
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Figure  23:  Dendrite  cell  size  in  ingots  1 and  2 (Al-4%  Cu) 
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Figure  26:  The  shape  of  the  mushy  zone  at  7 minutes  in  ingot 

(Al-4%  Cu). 
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Figure  27: 


Isotherm  detected  by  doping  ingot  4 


(Al-4%  Cu). 
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Figure  28:  Dendrite  cell  size  in  ingots  3,  4,  and  3 (Al-4%  Cu) . 
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The  macrosegregation  pattern  in  ingot  3 (Al-4%  Cu) . 


Figure  29: 
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Figure  30:  Shape  of  isotherms  in  ingots  3,  4,  and  5 (Al-4%  Cu) 

detected  by  doping. 
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Figure  31:  Uniform  composition  across  ingots  4 and  5 (Al-4%  Cu) . 
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Figure  32:  Position  of  the  liquidus  and  solidus  isotherms  in  the  Sn-15%  Pb  ingot. 
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Figure  33:  Secondary  dendrite  spacing  in  the  Sn-15%  Pb  ingot. 
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VI.  RESULTS  OF  COMPUTER  SIMULATIONS 


A.  Al-4%  Cu  Ingots 

Since  ingot  No.  1 exhibits  the  greatest  degree  of  macrosegregation, 
we  present  results  of  calculations  in  detail  for  that  ingot.  Figures  34 
and  35  are  plots  of  pressure  within  the  mushy  zone  calculated  assuming 
that  permeability  varies  with  fraction  liquid  according  to 

K = Yg^  (3) 

-7  2 

In  Figures  34,  y = 2 x 10  cm  and  pressure  increases  as  depth  down 
within  the  mushy  zone  increases.  If,  however,  permeability  within  the 

_ & 2 

mushy  zone  is  reduced  by  making  y = 8 x 10  cm  , flow  decreases  in  the 
lower  portion  of  the  mushy  zone  causing  a decrease  in  pressure  due  to 
friction  (Figure  35) . 

With  the  pressure  distribution  within  the  mushy  zone,  velocity  is 

calculated  according  to  Equation  (7).  It  is  of  more  interest,  however, 

to  examine  values  of  vVT/e  since  this  dimensionless  group  controls  the 

extent  of  macrosegregation  (see  Equation  (10)),  and  is  also  a measure  of 

channel  stability,  i.e.,  ''freckle1'  formation.  Figure  36  shows  that 

v-VT/e  is  greater  than  zero  so  there  should  be  no  channels  leading  to 

_ 6 2 

"freckles"  even  when  y is  made  as  great  as  10  cm  . In  the  central 
region  of  the  ingot  (Figure  36a)  the  effect  of  increasing  y is  to  reduce 
the  value  of  v*VT/e  and  to  increase  its  value  as  the  ingot  radius  is 
approached  (Figures  36b  and  c) . In  terms  of  macrosegregation,  this 
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effect  is  important  because  it  determines  whether  segregation  is  positive 
or  negative  at  ingot  center  or  ingot  radius. 

Figure  37  summarizes  results  of  calculations  for  values  of  y 

— 6 “8  2 7 2 

ranging  10  to  10  cm  . When  y < 5 x 10  cm  permeability  is  low  and 

"shrinkage  induced"  flow  predominates  over  "gravity  induced"  flow  leading 

to  negative  segregation  at  ingot  center  and  positive  segregation  at  ingot 
_ 6 2 

radius.  At  y = 10  cm  , permeability  is  high  and  "gravity  induced"  flow 

predominates  causing  the  dense  interdendritic  liquid,  rich  in  copper,  to 

flow  towards  the  ingot  center  thereby  leading  to  positive  segregation  at 

the  center  and  negative  segregation  at  ingot  radius. 

Agreement  between  calculations  and  experiment  for  ingot  1 appears 
_7  -7  2 

to  be  for  10  < y < 5 x 10  cm  . This  is  much  lower  than  the  value  of 

-9  -9  2 

3 x 10  to  9 x 10  cm  according  to  Figure  11,  but  secondary  arm 

spacing  is  also  significantly  less  (about  43  microns  versus  200  microns) . 

—8  2 

Flow  lines  are  shown  in  Figure  8 corresponding  to  y = 8 x 10  cm  . 

As  expected,  at  this  value  of  y,  "shrinkage  induced"  flow  (towards  the 

solidus  isotherm)  predominates,  and  there  is  no  flow  of  copper-rich 

liquid  towards  ingot  center;  hence,  segregation  is  negative  at  ingot 

center  and  positive  at  ingot  radius. 

In  the  central  regions  (r  < 2 cm.),  experimental  results  (Figure 

— 7 2 

37)  follow  the  calculated  curve  for  y = 5 x 10  cm  ; towards  the  ingot 

— 7 2 

radius  (2  £ r < 3.4  cm.),  there  is  better  agreement  with  y = 10  cm 
(Figure  37).  Thus,  it  appears  that  the  model  is  sufficiently  well  defined 
to  predict,  reasonably  well,  overall  variation  in  segregation  (i.e., 
whether  segregation  is  positive  or  negative  at  ingot  center  and  wall),  but 
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is  not  sensitive  enough  to  give  exact  agreement  with  experimental 
compositions  across  the  entire  ingot  cross  section. 

A thin  slice,  4 mm.  thick,  removed  from  across  the  ingot,  was 
radiographed  to  see  if  the  presence  of  porosity  could  account  for 
difference  between  experimental  results  and  calculations;  no  porosity  was 
detected  radiographically.  At  this  time,  it  appears  that  there  are  two 
reasons  for  the  discrepancy  between  calculated  and  experimental  results. 
Permeability  is  probably  not  isotropic  with  respect  to  flow  direction 
since  the  dendritic  structure  is  not  isotropic  nor  is  y uniform  throughout, 
because  dendrite  arm  spacing  is  not  constant  during  solidification  within 
the  mushy  zone.  The  second  reason,  and  probably  the  more  important  of  the 
two,  is  that  bulk  liquid  convection  influences  flow  within  the  mushy  zone 
so  that  the  boundary  condition  used  - Equation  (13)  - does  not  precisely 
describe  reality.  With  the  limited  knowledge  of  convection  within  liquid 
pools  of  ESR  ingots,  it  would  be  difficult  to  describe  this  boundary 
condition  precisely;  it  would  be  helpful,  however,  to  examine  the 
sensitivity  of  the  model  to  various  flow  patterns  imposed  at  the  liquidus 
isotherm  to  make  a judgment  on  the  usefulness  of  Equation  (13)  as  a 
boundary  condition. 

-7  2 

Using  y = 10  cm  , calculations  indicate  no  segregation  for 

ingot  2,  as  expected,  since  flow  lines  are  vertical  in  unidirectional 

solidification.  Calculations  are  also  compared  for  ingots  3,  4,  and  5 in 

-7  2 

Figure  38  assuming  y = 10  cm  . Again  the  model,  using  a single  value 
of  y,  does  not  predict  precisely  the  exact  segregation  pattern  (see 
Figures  29  and  31),  but  the  model  does  predict  less  segregation  in  ingots 


71 


4 and  5 than  in  ingot  3 which  is  expected  since  solidification  in  ingots 
4 and  5 more  nearly  approaches  unidirectional  solidification  than  does 
solidification  in  ingot  3 (Figure  30). 


B.  Sn-15%  Pb  Ingot 

As  in  the  case  of  Al-4%  Cu  ingots,  agreement  between  experimental 

and  calculated  results  is  not  precise,  but  the  computer  model  does 

predict  the  overall  segregation  pattern,  across  an  ingot,  reasonably 

—6  2 

well.  Figure  39  indicates  that  calculations  using  Y = 10  cm  agree 

-7  2 

more  closely  with  experimental  results  than  when  y - 5 x 10  cm  . On 

this  basis  alone,  it  appears  that  the  permeability  (at  a given  fraction 

liquid)  is  greater  in  Sn-15%  Pb  alloys  than  in  Al-4%  Cu  alloys  by  about 

one  order  of  magnitude.  Measured  values  of  permeability  in  lead  rich 

1 2 

Pb-Sn  alloys  are  reported  to  be  10  - 10  greater  than  the  value  of 

permeability  reported  for  aluminum  alloys . Assuming  similar 

behavior  of  tin-rich  and  lead-rich  alloys,  it  is  reasonable,  therefore, 

—6  2 — 7 2 
to  use  y = 10  cm  in  the  Sn-15%  Pb  ingot  contrasted  to  y = 10  cm 

in  the  Al-4%  Cu  ingot. 

In  the  Sn-15%  Pb  ingot,  flow  is  extensive  as  indicated  by  the  values 
of  vVT/e  shown  in  Figure  40.  In  most  of  the  central  region  (Figure  40a), 
vVT/e  is  negative  and  drops  below  -1  for  r < 0.8  cm.  Away  from  the 
center  (Figures  40b  and  c) , vVT/e  is  sometimes  negative  but  predominantly 
positive;  it  never  drops  below  a value  of  -1.  Figure  9 shows  that  in  much 
of  the  ingot,  interdendritic  flow  is  from  colder  to  hotter  isotherms. 
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Obviously,  in  this  ingot  " gravity  induced"  flow  predominates  and 
consequently  segregation  is  not  only  positive  towards  the  center,  but 
there  is  a high  degree  of  channeling  and  copious  freckling.  To  sum  up, 
the  analytical  model  predicts  the  occurrence  of  "freckles"  which  were  also 
found  in  the  experimental  ingot  of  Sn-15%  Pb  alloy  (Figure  7). 

Work  is  currently  underway  to  predict  conditions  necessary  to  avoid 
formation  of  "freckles"  in  ESR  ingots,  using  the  model  alloy  and  the 
computer  program  together. 


PRESSURE  (P- 


DISTANCE  UP  FROM  SOLIDUS  (z-zs),  cm 

- 7 2 

Figure  34:  Pressure  in  ingot  1 calculated  with  y = 2 x 10  cm  , 


(T*  VT)/€ 


COMPOSITION,  % Cu 


COMPOSITION,  % Cu 
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Figure  38:  Calculations  of  segregation  in  ingots  3,  4,  and  5 with 

Y = 10"  7 cm2. 
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qd%  ‘NOIllSOdlNOO 


RADIUS,  cm 

Figure  39:  Calculation  of  segregation  in  the  Sn-15%  Pb  ingot. 
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VII.  CONCLUSIONS 


1.  Severe  macrosegregation  in  ESR  ingots,  such  as  "freckles"  and 
other  "channel"- type  segregates,  as  well  as  overall  surface  to 
center  variations  in  composition,  can  be  produced  and  studied 
using  laboratory-scale  apparatus. 

2.  Calculations  of  macrosegregation,  using  a theoretical  model,  agree 
well  with  experimental  results  considering  the  number  of  variables 
involved.  The  model  calculates  surface  to  center  variations  in 
composition  using  as  input  the  temperature  distribution  within  the 
"mushy  zone"  of  a solidifying  ESR  ingot.  The  model  also  predicts 
the  formation  of  "freckles"  and  "channel"- type  segregates. 

3.  Heat  flow  conditions  which  tend  to  make  deeper  isotherms  within 
the  "mushy  zone"  accentuate  segregation  in  ESR  ingots.  On  the 
other  hand,  in  ingots  with  low  cooling  rates,  permeability  within 
the  partially  solidified  alloy  is  high  (because  dendrite  arms  are 
more  widely  spaced) , and  interdendritic  flow  of  solute-rich 
liquid  is  extensive,  leading  to  severe  macrosegregation. 
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APPENDIX  A - ANALYSIS  OF  FLOW  OF  INTERDENDRITIC  LIQUID 

In  analyzing  the  flow  within  the  mushy  zone,  the  volume  element 
taken  is  large  enough  so  that  fraction  solid  within  it  is  exactly  the 
local  average,  but  small  enough  so  that  it  can  be  treated  as  a 
differential  element.  This  is  the  typical  differential  volume  element 
chosen  in  two-phase  flow.  For  solidification  problems  it  contains 
several  dendrite  arms.  To  analyze  macrosegregation  many  of  the  details 
of  solute  redistribution  on  the  microscale  must  be  described.  For 
dendritic  freezing  such  as  the  model  alloys  selected  for  this  study,  we 
can  assume  the  following  : 

1.  Liquid  composition  and  temperature  within  the  differential 
are  uniform. 

2.  Equilibrium  exists  at  the  solid-liquid  interface  so  that  the 

equilibrium  partition  ratio,  k,  applies  and  k = C*/CT  where  C*  and  C 

s L s L 

are  the  compositions  of  the  interface  solid  and  liquid  phases, 
respectively,  and  their  values  are  from  the  equilibrium  phase  diagram 
at  the  temperature  of  the  differential  element. 

3.  Diffusion  in  the  solid  is  negligible. 

4.  Only  mass  flow  in  or  out  of  the  differential  element  by 
convection  is  considered.  Mass  flow  in  or  out  of  the  differential 
element  by  diffusion  in  the  liquid  or  mass  input  by  flow  of  settling 
solid  phase  is  neglected. 
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5.  Density  of  solid,  p , during  solidification  is  constant.  This 

s 

is  a very  close  approximation  for  Al-4%  Cu  and  Sn-15%  Pb  systems  used  as 
model  alloys  herein. 

6.  No  pore  formation;  i.e.,  g + g = 1 where  g and  g are  volume 

S t S Li 

fraction  solid  and  liquid,  respectively. 

Conservation  of  mass  within  the  differential  element  is 


at 


(ps8s  + PLgL) 


■v-p 


8t^ 


L L 


(20) 


where  t is  time  and  v is  velocity  of  the  interdendritic  liquid.  With  g 

o 

plus  g =1  and  p constant,  Equation  (20)  becomes 
L s 


^pl  “ ps)  at-  + 8l  at-  = -v'pl8lv 


(21) 


In  the  mushy  zone  we  apply  D'Arcy's  Law  for  flow  through  porous 
media.  Therefore 

V = - -§-  (VP  + p g)  (1) 

where  K is  the  permeability*,  P is  pressure,  and  y is  the  liquid  viscosity. 
A relationship  of  the  form 

K = yg l (3) 


* K is  actually  the  "specific  permeability"  as  defined  in  Geiger  and 
(39) 


Poirier 


but  is  called  permeability  in  this  report 
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has  been  used  by  Mehrabian  et  al^^  to  study  macrosegregation  based  on 
the  experimental  work  of  Flemings  and  coworkers 9 . 

By  substituting  Equations  (1)  and  (3)  into  Equation  (21),  we  get 


(Pi 


P ) 


9§i 


3Pi 


s'  at 


+ g 


l at 


2 

YgLPL 

( — — VP  + 


2 2 

Y8lpl  -> 


g ) 


(22) 


During  dendritic  solidification  the  interdendritic  liquid  is  a 

function  of  temperature  only,  i.e.,  C versus  T is  given  by  the  phase 

-L 

diagram.  This  follows  from  assumptions  (1)  and  (2).  Since  C is  a 

Li 

function  of  T only,  we  can  write 


)( 


3t 

at 


e 

m 


(23) 


where  m is  the  slope  of  the  liquidus  from  the  equilibrium  phase  diagram, 

and  £ is  the  local  cooling  rate  of  the  differential  element. 

Based  upon  the  assumptions  listed  above,  the  "local  solute 

(13) 

redistribution  equation"  has  been  derived v . It  gives  volume  fraction 
liquid  as  a function  of  the  local  composition  of  the  interdendritic 
liquid  (and  hence,  as  a function  of  temperature  during  solidification); 
it  is 


)(  1 + 


(9) 


From  the  chain  rule,  we  have 


3C, 


)( 


5 

at 


(24) 
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and  by  substituting  Equations  (23)  and  (9)  into  Equation  (24),  we  get 


9sl 

3t 


k)(1  + 


v-VT,  £l  t 

e CL  m 


(25) 


The  density  of  the  interdendritic  liquid,  p , is  a function  of 

Li 

concentration  and  temperature.  But  since  C = f(T),  then  p is  a function 

Li  JL 

of  CL  only  (Figures  13b  and  14b).  Therefore, 


9pl 

3t 


dp, 

dC, 


)( 


9C, 

at 


(26) 


where  (dp  /dC  ) is  the  slope  of  a plot  of  p versus  C . 

Li  Li  Li  L 

By  substituting  Equations  (25)  and  (26)  into  Equation  (22),  and 
applying  D'Arcy’s  Law,  we  get 
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( — — VP  + 
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2 2 
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Yg, 

ye 


(VP  + pLg) 


gT  p 

VT  ] 

Li 
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(27) 


Equation  (27)  can  be  solved  along  with  appropriate  boundary 
conditions  to  give  pressure  within  the  mushy  zone  of  a solidifying  ingot. 
For  application  to  our  analytical  model  - specifically  for  cylindrical 
ESR  ingots  - we  expand  Equation  (27)  into  cylindrical  coordinates  (z,r), 


and  use 
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ii|  + AM  + A + b»p  + c , 0 _ 
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where 
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and 


a = 
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(1  - k)C, 


(10) 
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APPENDIX  B - METHOD  OF  SOLUTION 


1.  Finite  Difference  Approximations 

Equation  (10),  along  with  the  boundary  conditions.  Equations  (11) 
(14),  are  solved  using  a computer  program  based  on  finite  difference 
approximations  of  the  derivatives. 


Interior  Nodes 


Figure  41a  depicts  an  interior  node  within  the  mushy  zone  at  a 
pressure  PQ  surrounded  by  four  nodes  at  pressures  P^,  P^,  P^,  and  P^ . 
Node  spacing  in  the  r-direction  is  Ar  = h and  in  the  z-direction  is 
Az  = k.  The  finite  difference  approximations  for  the  pressure 
derivatives  are 


8P 

3r 


9 


similarly 


9p  h h 

3z  - 2k 


also 
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8r 
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'■3r; 


P2  + 2Pq  + P4 


and 
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(28) 


(29) 


(30) 


(31) 
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The  coefficients  A,  B,  and  C are  evaluated  by  assigning  a temperature 
at  each  node,  thus  determining  C , p , a,  y,  £,  and  (dp  /dC  ) . To  initiate 

J-j  ii  L L 

calculations,  values  of  g^  at  each  node  are  approximated  by  the  Scheil 
equation  (i.e.,  Equation  (9)  with  v = 0 and  3=0);  by  an  iteration 
technique,  values  of  g are  "updated"  as  calculations  proceed. 

ij 

Boundary  Nodes  Along  the  Liquidus 

The  liquidus  isotherm  does  not  necessarily  pass  through  the 
location  of  nodal  points  so  "boundary  nodes"  are  defined  by  determining 
the  node  (above  or  below  the  isotherm)  which  is  closer  to  the  isotherm. 

That  node  is  then  considered  to  be  a "boundary  node".  Specifically  in 
Figure  41b,  the  boundary  nodes  are  2,  3,  and  6 because  they  are  closer  to 
the  liquidus  isotherm  than  nodes  1,  4,  and  5,  respectively. 

The  assumption  of  no  convection  in  the  liquid  pool  is  made  so  that 
the  pressure  along  the  liquidus  is  given  by  the  metallostatic  head  above 
the  liquidus.  Thus,  taking  node  3 as  example, 

P3  ' Pa+PL0St  <33> 

where  P is  the  pressure  at  the  top  of  the  liquid  pool  (1  atmosphere  is 
a 

used  in  our  calculations),  i is  the  height  of  liquid  metal  above  node  3 
and  p^Q  is  the  density  of  the  liquid  metal  pool. 

If  the  liquidus  isotherm  is  steep  enough,  then  locating  "boundary 
nodes"  is  done  by  deciding  which  of  two  nodes  at  a given  z-value  is  the 
closer  rather  than  deciding  which  of  two  nodes  at  a given  r-value  is 


92 


closer  as  depicted  above.  Obviously  as  the  grid  spacing  decreases,  then 
the  positions  of  the  "boundary  nodes"  more  nearly  approximate  the  liquidus 
isotherm.  For  ingot  1,  calculations  were  done  with  k = 0.1  cm  and  h = 0.2 
cm;  in  ingot  2,  k = 0.2  cm  and  h = 0.2  cm.  In  ingots  3,  4,  and  5 
k = 0.03  cm  and  h = 0.2  cm.  Calculations  using  k = 0.14  cm  and  h = 0.25 
cm  were  made  for  the  Sn-15%  Pb  ingot. 


Boundary  Nodes  Along  the  Solidus 


At  the  solidus  isotherm,  alloy  solidification  is  completed  by 
eutectic  liquid  solidifying  to  form  eutectic  solid  of  the  same 
composition.  If  no  pores  form,  then  the  flow  of  interdendritic  liquid 
must  equal  the  rate  at  which  the  solidifying  eutectic  liquid  is 
shrinking.  Hence, 


v 


SE 


- P 


U, 


LE 


(14) 


where  Ug  is  the  eutectic  isotherm  velocity  normal  to  the  isotherm.  The 
components  of  Equation  (14)  are 
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Equations  (33)  and  (34)  are  in  terms  of  velocity;  it  is  necessary  to 

transform  them  into  statements  which  describe  pressure  gradient  along  the 

boundary.  This  is  done  by  using  the  components  of  D’Arcy’s  equation, 

Equation  (1),  with  p = p and  gT  = gTt?,  as  follows: 
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By  combining  Equations  (33)  and  (35)  and  Equations  (34)  and  (36),  we 
express  the  components  of  the  pressure  gradient  at  the  solidus  isotherm. 
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(38) 


To  formulate  these  pressure  gradients  into  the  finite-difference 
computer  program,  we  must  first  locate  the  ’’boundary  nodes"  along  the 
solidus  in  a manner  similar  to  locating  "boundary  nodes"  for  the  liquidus. 
Figure  41c  shows  a group  of  nodes  in  the  vicinity  of  the  solidus  isotherm. 
The  "boundary  nodes”  are  selected  by  deciding  if  is  closer  to  the 
isotherm  than  is  P^,  P^  closer  than  P^,  and  P^  closer  than  P^ • In  this 
example,  P^,  P , and  P^  are  selected  as  "boundary  nodes".  Then  the 
finite  difference  approximations  for  Equation  (10)  are  applied  to  each 
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node  directly  above  the  "boundary  nodes".  Node  Pq  is  such  a node.  The 
derivatives  appropriately  take  into  account  the  components  of  pressure 
gradient  at  the  solidus  (Equations  (37)  and  (38)).  These  derivatives, 
in  finite  difference  form,  the 
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This  is  all  done  in  the  same  manner  as  for  an  interior  nodes,  except  now 
we  use  Equation  (38)  and  write 
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(43) 


Equation  (43)  is  substituted  into  Equations  (41)  and  (42)  to  obtain  the 
z-derivates  for  node  "0".  In  assembling  the  program,  therefore,  there  is 
an  equation  for  each  of  the  nodes  just  above  each  "boundary  nodes". 

There  are  no  equations  for  the  "boundary  nodes",  per  se. 

If  the  solidus  isotherm  is  steep,  then  "boundary  nodes"  are 
selected  by  choosing  the  closer  node  to  the  left  or  right  of  the  isotherm. 
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Then  derivatives  in  the  r-direction  are  modified  and  Equation  (37)  is 
invoked  in  a manner  similar  to  that  described  above  for  Equation  (38) . 

Boundary  Nodes  at  the  Centerline  and  Ingot  Radius 

By  symmetry,  v^  = 0 or  3P/3r  = 0 at  the  centerline.  At  ingot  radius, 
vf  = 0 so  along  this  boundary  9P/3r  = 0 also.  These  boundary  conditions 
are  incorporated  into  the  computer  program  by  merely  setting  pressure  at 
the  centerline  nodes  equal  to  the  pressure  of  the  nodes  directly 
connected  to  the  centerline  nodes.  Similarly,  the  pressure  at  a node 
along  the  ingot  radius  is  set  equal  to  its  connecting  node. 

2.  Computer  Program 

When  approximating  equations  are  written  for  all  of  the  interior 
and  boundary  nodes,  there  results  a set  of  n equations  which  must  be 
solved  for  n unknown  values  of  pressure.  Since  n is  in  the  range  of 
500  - 1000,  the  set  of  equations  are  not  solved  directly,  but  rather  an 
iteration  scheme  is  used.  Figure  15  is  a flow  chart  of  the  computer  code. 
First  pressure  at  all  nodes  is  approximated  by  merely  calculating  pressure 
due  to  metallostatic  pressure.  Then,  node  by  node,  the  computer  code  is 
directed  to  solve  for  nodal  pressures  using  the  "first  guess"  pressures 
at  connecting  nodes;  then  "second  guess"  values  of  pressure  are 
calculated  and  assigned  to  every  node.  Then  the  entire  process  is 
repeated  200  times  (value  of  MAXIT) . 
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For  the  first  200  iterations,  the  volume  fraction  liquid  at  each 
node  is  that  calculated  by  the  Scheil  equation  corresponding  to  the  nodal 
temperature.  In  Figure  15  this  brings  us  to  step  7. 

With  pressure  distribution  in  the  mushy  zone  calculated,  velocity 
components  at  all  nodes  can  then  be  calculated  by  using  DfArcyTs  Law  - 
Equation  (1)  - which  relates  velocity  to  pressure  gradient.  The  "local 
solute  redistribution  equation"  - Equation  (9)  - is  then  integrated, 
node  by  node,  down  a given  column  from  the  liquidus  to  the  solidus  in 
order  to  obtain  updated  values  of  Equation  (15)  is  solved  to  obtain 

local  average  composition.  With  these  updated  values  of  g (step  8 in 

Lt 

Figure  15),  pressure  is  recalculated  for  200  more  iterations  (MAXIT=400) , 
bringing  the  total  iterations  up  to  400  when  a new  set  of  g values  are 
obtained.  The  process  is  repeated  for  another  200  iterations. 

A total  of  600  iterations  was  selected  because  a trial  run 
showed  that  from  500  to  700  iterations  only  a negligible  change  in  local 
average  compositions  was  noted  when  the  computer  code  was  used  to 
simulate  macro segregation  in  ingot  1. 
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Figure  41:  Nodes  used  in  the  finite  difference  approximations  of 

pressure  within  the  mushy  zone:  (a)  interior  nodes; 

(b)  the  nodes  at  the  liquidus  isotherm;  (c)  nodes  at 
the  solidus  isotherm. 
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first,  employed  for  Al-4%Cu,  is  a small  scale  ESR  unit  which  consists  of 
a DC  power  source,  a water  cooled  mold,  consumable  electrodes,  and  slag 
layer  as  in  conventional  ESR  processing.  The  second  type  of  apparatus 
is  one  which  simulates  the  solidification  of  an  ESR  ingot  but  does  not 
employ  slag.  The  unit,  used  for  Sn-15%Pb,  comprises  a source  of  melt 
droplets,  a cooled  mold,  and  a heat  source  to  simulate  the  heat  input 
of  the  ESR  process. 

Results  of  macrosegregation  in  a series  of  five  Al-4%Cu  experimental 
ingots  show  that  severe-localized  segregates,  as  sometimes  found  in 
large  commercial  ingots,  cannot  be  produced  using  this  alloy  cast  in  the 
small-scale  ESR  molds  (7.6  cm  in  diameter).  Surface  to  center 
compositional  variations  are  produced  enabling  a comparison  with 
calculations  using  macrosegregation  theory. 

In  order  to  study  a wider  range  of  segregation  problems  encountered  in 
ESR  ingots,  recent  efforts  are  directed  toward  solidifying  Pb-15%Sn 
alloy  in  the  simulated  ESR  apparatus.  Initial  results  obtained  are 
very  encouraging.  Severe  macrosegregation  is  produced  in  this 
apparatus.  In  an  ingot  (8.2  cm  diameter),  the  composition  varies  from 
about  26%  Pb  in  the  center  to  a minimum  of  6%Pb  slightly  beyond 
midradius.  Additionally,  there  are  ’’freckles"  at  approximately 
midradius  and  an  abnormally  high  fraction  of  eutectic  constituent 
(similar  to  "V-segregates")  in  the  center. 

The  analytical  work  is  based  upon  an  analysis  of  the  slow  creeping 
flow  of  liquid  between  dendrite  arms.  The  flow  is  driven  by  (1) 
gravity  induced  convection,  and  (2)  solidification  shrinkage.  Agreement 
between  calculations  of  macrosegregation  and  the  macrosegregation  in  the 
experimental  ingots  is  excellent  considering  the  large  number  of 
variables  involved.  The  theory  correctly  predicts  the  presence  of 
"freckles'*  found  in  the  Sn-15%Pb  ingot. 
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